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ABSTRACT 

Context. Recently, evidence for the presence of weak magnetic fields in Herbig Ae/Be stars has been found in several 
studies. 

Aims. We seek to expand the sample of intermediate-mass pre-main sequence stars with circular polarization data to 
measure their magnetic fields, and to determine whether magnetic field properties in these stars are correlated with 
mass-accretion rate, disk inclination, companions, silicates, PAHs, or show a correlation with age and X-ray emission 
as expected for the decay of a remnant dynamo. 

Methods. Spectropolarimetric observations of 21 Herbig Ae/Be stars and six debris disk stars have been obtained at the 
European Southern Observatory with FORS 1 mounted on the 8 m Kueyen telescope of the VLT. With the GRISM 600B 
in the wavelength range 3250-6215 A we were able to cover all hydrogen Balmer lines from H/3 to the Balmer jump. In 
all observations a slit width of O'.'l was used to obtain a spectral resolving power of i? ~ 2000. 

Results. Among the 21 Herbig Ae/Be stars studied, new detections of a magnetic field were achieved in six stars. 
For three Herbig Ae/Be stars, we confirm previous magnetic field detections. The largest longitudinal magnetic field, 
(Bz) = — 454±42 G, was detected in the Herbig Ae/Be star HD 101412 using hydrogen lines. No field detection at a 
significance level of 3a was achieved in stars with debris disks. Our study does not indicate any correlation of the 
strength of the longitudinal magnetic field with disk orientation, disk geometry, or the presence of a companion. We 
also do not see any simple dependence on the mass-accretion rate. However, it is likely that the range of observed field 
values qualitatively supports the expectations from magnet ospheric accretion models giving support for dipole-like field 
geometries. Both the magnetic field strength and the X-ray emission show hints of a decline with age in the range of ~2- 
14Myrs probed by our sample, supporting a dynamo mechanism that decays with age. However, our study of rotation 
does not show any obvious trend of the strength of the longitudinal magnetic field with rotation period. Furthermore, 
the stars seem to obey the universal power-law relation between magnetic flux and X-ray luminosity established for the 
Sun and main-sequence active dwarf stars. 

Key words, polarization - Stars: pre-main-sequence - Stars: circumstellar matter - Stars: magnetic fields - Xrays: stars 
- Stars: coronae - Stars: activity 



1. Introduction 

Magnetic fields are important ingredients of the star forma- 
tion process (McKee & Ostriker l2007p . Models of magnet- 
ically driven accretion and outflows (e.g., Shu et al. 120001 
Shu et al. I1995P successfully reproduce many observational 
properties of low-mass pre-main sequence stars, the clas- 
sical TTauri stars (cTTS). Indirect observational evidence 
for the presence of magnetic fields in these stars is seen 



* Based on observations obtained at the European Southern 
Observatory, Paranal, Chile (ESO programmes 077.C-0521(A) 
and 081.C-0410(A)). 



in strong X-ray, FUV, and UV emission (e.g., Feigelson & 
Montmerle 119991 Brown et al. I1985p . The high-energy ra- 
diation of young stars related to magnetic processes may 
be critical for the annealing and melting of some refractory 
grain species (Shu et al. l200ip . for the heating intermediate- 
depth portions of disks around young stars fNaiita r2004p 
and for mass loss from planets forming in those disks 
(Penz et al. I2008p . The irradiation- induced heating medi- 
ated by the magnetic activity enlarges the zone in circum- 
stellar disks where chemistry may proceed to assemble pre- 
biotic materials which could produce life-bearing worlds. 
Therefore, understanding the interaction between central 
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stars, their magnetic fields and protoplanctary disks is cru- 
cial to reconstruct the Solar System's history, and to ac- 
count for the diversity of exo-planetary systems. 

Direct magnetic field measurements of cTTS (e.g. 
Johns-KruU I2007P corroborate the scenario outlined above 
for the formation of low-mass stars, but the picture is 
less clear for higher-mass stars. The presence of protoplan- 
etary disks around intermediate-mass pre-main sequence 
(PMS) stars, the Her big A e/Be stars (e.g. Herbig UMUi 
Finkenzeller & Mundt fTM4l The et al fTOM]) . has been es- 
tablished through their thermal emission in the IR, and 
more recently by direct imaging with polarimetric and in- 
terferometric techniques (e.g., Perrin et al. 120041 Monnier 
et al. l2005p . Dusty disk models explain the IR excess emis- 
sion (e.g. DuUemond et al. I200ip and magnetospheric ac- 
cretion models describe the line profiles of Herbig Ac/Be 
stars (MuzeroUe et al. I2004p . in close analogy to cTTS. 
However, the disk lifetime appears to be shorter for Herbig 
Ae/Be stars (Uzpen et al. 120091 Carpenter et al. l2005p and 
within the Herbig Ae/Be class there is a trend of more 
massive stars dispersing their disks more rapidly (Alonso- 
Albi et al. I2009p . Herbig Ae stars are closer analogs to 
TTS than Herbig Be stars. Bipolar outflows or associated 
Herbig-Haro knots are now known for six optically visible 
Herbig Ae stars drawn from coronographic imaging surveys 
(e.g. Melnikov et al.[2008j). Due to the probable role of mag- 
netic fields in launching and coUimating jets, this implies 
that a significant fraction of Herbig Ae stars should have 
measurable magnetic fields, which could persist through 
much, if not all, of the star's PMS lifetime. X-ray activ- 
ity is known for a number of Herbig Ae stars (Hamaguchi 
et al. [20051 Feige lson et al. [20031 Swart z et al. [5005l Stelzer 
et aL l^OMI^DTOl Skinner et al. [^D(H|) . suggesting the pres- 
ence of magnetic fields. Recent advances in instrumenta- 
tion have resulted in the first magnetic field measurements 
for Herbig Ae stars, derived from spectropolarimetry. A 
few stars have magnetic field strengths derived from cir- 
cular polarization measurements of a few hundred Gauss 
or near 100 G (Hubrig et al. Wade et al. 20051 IW71 

Hubrig et al. I2006al I2007ap . Other objects seem to have 
either smaller average field strengths, or exhibit significant 
variability. 

In our previous studies, we reported detections at a level 
higher than 3a for three out of seven Herbig Ae stars ob- 
served with FORS 1 (Hubrig et al. I2004al I2006al I2007ap . 
The results of the spectropolarimetric observations of the 
two Herbig Ae stars HD 31648 and HD 190073 with FORS 1 
at a resolution of i? w 4000 turned out to be especially re- 
markable (Hubrig et al. l2006a[ Hubrig et al. l2007ap . These 
two stars possess almost the same effective temperature, 
but are different in luminosity and projected rotation ve- 
locity. In the Stokes V spectra of HD 190073, the profiles ex- 
hibit a number of blucshiftcd local absorption components, 
whereas in HD 31648, only one blucshiftcd and one red- 
shifted feature are observed in both H&K lines of the Ca ii 
doublet. It is noteworthy that recent Goddard Fabry-Perot 
narrow band images revealed the counterjet and bright HH 
knots in HD 31648 (Stecklum et al. 2009, in preparation). 
This star demonstrates notable emission in the H/3, H7, 
and H(5 lines, which indicates the presence of a significant 
stellar wind. The profiles of these lines arc of P Cyg-type 
with very deep blucshiftcd absorptions. 

In this work we seek to expand the sample of 
intermediate-mass PMS stars with circular polarization 



data used to derive stellar magnetic fields. We then investi- 
gate whether the magnetic field properties of these stars are 
correlated with mass-accretion rate, companions, silicates, 
or PAHs. We also search for the first time for a correlation 
between field strength and age or X-ray emission, such as 
expected for a decaying fossil field. Further, we examine 
the inclination dependence of the polarization data, and 
explore the dependence of the strength of the longitudinal 
magnetic field on stellar rotation. 

2. Observations and data reduction 

The observations reported here were carried out on May 22 
and 23, 2008 in visitor mode at the European Southern 
Observatory with FORS 1 mounted on the 8 m Kueyen 
telescope of the VLT. This multi-mode instrument is 
equipped with polarisation analyzing optics comprising 
super-achromatic half-wave and quarter-wave phase re- 
tardcr plates, and a WoUaston prism with a beam diver- 
gence of 22" in standard resolution mode. In 2007, a new 
mosaic detector with blue optimized E2V chips was in- 
stalled in FORSl. It has a pixel size of 15 fim (compared 
to 24 fim for the previous Tektronix chip) and a higher ef- 
ficiency in the wavelength range below 6000 A. To achieve 
the highest possible signal-to- noise (S/N) ratio ~ as is re- 
quired for accurate measurements of stellar magnetic fields 
- we used the (200kHz, low, 1x1) readout mode, which 
makes it possible to achieve a S/N ratio of about 1000- 
1200 with only one single exposure. With the GRISM600B 
in the wavelength range 3250-6215 A we were able to cover 
all hydrogen Balmer lines from H/? to the Balmer jump. 
In all observations, a slit width of 0'.'4 was used to obtain 
a spectral resolving power of « 2000. Usually, we took 
three to five continuous series of two exposures for each 
object, leading to signal-to-noise ratios of typically 3000 to 
4000. For the faintest Herbig Ae star in our sample, VV Ser, 
only two series were taken due to the rather long exposure 
time for each sub-exposure. More details on the observ- 
ing technique with FORSl can be found elsewhere (e.g., 
Hubrig et al. I2004al I2004c[ and references therein). 

The mean longitudinal magnetic field is the average over 
the stellar hemisphere visible at the time of observation of 
the component of the magnetic field parallel to the line of 
sight, weighted by the local emergent spectral line intensity. 
Its determination is based on the use of the equation 
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.gofreA^ 1 dl 



ATmieC^ I dX 



(1) 



where V is the Stokes parameter which measures the 
circular polarization, / is the intensity observed in unpo- 
larized light, (70ft is the effective Lande factor, e is the elec- 
tron charge, A is the wavelength, rrie the electron mass, c 
the speed of light, dl/dX is the derivative of Stokes /, and 
(Bz) is the mean longitudinal magnetic field. To minimize 
the cross-talk effect, we executed the sequence -t-45— 45, 
-1-45—45, -1-45—45 etc. and calculated the values V/I using: 
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r + r 



(2) 



Q=+45° 



where a denotes the position angle of the retarder wave- 
plate and f° and /° are ordinary and extraordinary beams. 
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respectively. Stokes / values were obtained from the sum of 
the ordinary and extraordinary beams, which are recorded 
simultaneously by the detector. To derive (Bz), a least- 
squares technique was used to minimize the expression 

2 >^ {y^ - (B,) X, - hf 

X =2^ -2 (3) 

i ^ 

where, for each spectral point i, yi — {V/I)i, Xi = 

~^^ri? (V-^ ^ d//dA)i, and 6 is a constant term that, 
assuming that Eq. [l]is correct, approximates the fraction 
of instrumental polarization not removed after the applica- 
tion of Eq.[2]to the observations. During the commissioning 
of FORS 1, this instrumental polarization term was found 
to be wavelength independent. We assume that the only 
source of uncertainty in our field measurements is from the 
photon count statistics of the observations. Therefore, the 
longitudinal field uncertainty is obtained from the formal 
uncertainty of the linear regression. For each spectral point 
i, the derivative of Stokes / with respect to the wavelength 
was evaluated following 

(f\ . (4) 

where Ni is the photon count at wavelength A^. Since noise 
strongly influences the derivative, we interpolate the data 
after spectrum extraction with splines. Our experience from 
a study of a large sample of magnetic and non-magnetic 
Ap and Bp stars (Hubrig et al. I2006b|) revealed that this 
regression technique is very robust and that detections with 
(5z)>3(T result only for stars possessing magnetic fields. 

The list of the studied stars is presented in Table [TJ In 
the four columns we give the object name, another identifier 
from SIMBAD, the visual magnitude and the spectral type. 
The star HD 47839 of spectral type 07 V was included in 
our sample since it is classified as a pre-main sequence star 
in the SIMBAD database, probably due to its close proxim- 
ity to the Cone Nebula. Markova et al. I|2004p consider this 
star as a Galactic 0-type star with a mass of 32 M0 and 
Teff=37500K. In Table [5] we Hst various properties of the 
studied stars compiled from the literature, such as wsini, 
mass-accretion rate, features in the mid-IR spectra, the 
presence of a companion, effective temperature, luminosity, 
and age. The mass was deduced from the position of the star 
in the H-R diagram (see Sect. [S]). The mass- accretion rates 
were determined by Garcia Lopez et al. (|2006p from the lu- 
minosity of the Br7 line seen in emission in Herbig Ae/Be 
stars. Mid-IR features such as crystalline silicate grains in- 
dicate thermal dust processing subsequent to the formation 
of the star-disk system, while PAHs (polycyclic aromatic 
hydrocarbons) are abundant in interstellar and eircumstel- 
lar environments and can be important constituents in the 
energy balance of those environments as sources of photo- 
electrons that heat the gas component (Kamp & DuUemond 
I2004p . Since the existence of a link between the presence 
of magnetic fields and the disk characteristics remains un- 
proven, we decided to include in this study also the infor- 
mation on mid-IR features. Additionally, in Tables [T] and 
[5] we added five Herbig Ae/Be stars for which magnetic 
fields were detected in previous studies by various authors. 
The magnetic field for HD 31648 was detected by Hubrig 
et al. (|2006bp . while V380Ori and BFOri were studied by 



Table 1. Target stars for which spectropolarimetric data 
were obtained during our observing run. 



Object 


Other 


V 


Spectral 


name 


identifier 




type 


PDS2 


CPD-53 295 


10.7 


F2 


HD 47839 


15Mon 


4.7 


07 Ve 


HD 95881 


Hen 3-554 


8.3 


Al HI 


HD 97048 


CUCha 


8.5 


AOpshe 


HD 97300 


CED 112 IRS 3 


9.0 


B9 V 


HD 100453 


CD-53 4102 


7.8 


A9 Ve 


HD 100546 


KRMus 


6.7 


B9 Vne 


HD 101412 


V1052 Cen 


9.0 


B9.5 V 


HD 135344B 


SAO 206462 


7.8 


F4-F8 


HD 139614 


Hen 3-1086 


8.3 


A7 Ve 


HD 144432 


Hen 3-1141 


8.2 


A9/F0 V 


HD 144668 


HR5999 


7.0 


A7 IVe 


HD 150193 


V2307 Oph 


8.9 


Al Ve 


HD 152404 


AK Sco 


9.1 


F5 Ve 


HD 158643 


51 Oph 


4.8 


AO V 


HD 163296 


MWC 275 


6.9 


Al Ve 


HD 169142 


MWC 925 


8.2 


B9 Ve 


VVSer 


FMC39 


11.6 


A2e 


HD 176386 


CD-37 13023 


7.3 


B9 IVe 


HD 179218 


BD-f 15 1991 


7.2 


B9e 


HD 190073 


MWC 325 


7.8 


A2 IVpe 


Stars with debris disks 


HD 9672 


49 Get 


5.6 


Al V 


HD 39060 




3.9 


A6V 


HD 109573 


HR4796 


5.8 


AO V 


HD 164249 


CD-51 11312 


7.3 


F5 V 


HD 172555 


HR7012 


4.8 


A7 V 


HD 181327 


CD-54 8270 


7.0 


F5/F6V 


Magnetic Herbig Ae/Be stars from previous studies 


HD 31648 


MWC 480 


7.7 


A3pshe 


HD 104237 


Hen 3-741 


6.6 


A:pe 


HD 200775 


MWC 361 


7.4 


B2 Ve 


V380 Ori 


MWC 765 


10.7 


AO 


BF Ori 


BD-061259 


10.3 


A5II/IIIev 


Standard Ap star with a weak 


mag 


netic field 


HD 162725 


V951 Sco 


6.4 


Ap 



Notes: Spectral types and visual magnitudes are taken from 
SIMBAD. 



Wade et al. (|2005p . The detection of a weak magnetic field 
in HD 104237 was reported by Donati et al. (|T997|) and 
the study of the magnetic field of HD 200775 was carried 
out by Alecian et al. (|2008p. W e note, however, that, as was 
shown by Hubrig et al. (j2007bp . the detected magnetic field 
of HD 31648 is predominantly circumstellar (CS) and that 
the strength of the photospheric magnetic field remains un- 
known. Also, the field detection in HD 104237 of the order 
of ^^50 G was just marginal, and was not confirmed in the 
follow-up study by Wade et al. (|2007p . 

Since caution is called for in the detection of mag- 
netic fields with low resolution spectropolarimeters, from 
time to time we observed magnetic Ap/Bp stars with well- 
known variation curves. These observations confirm that 
our measurements are usually in good agreement with mea- 
surements obtained with other spectropolarimeters (see 
e.g. Hubrig et al. I2004bp . During our observing run, we 
also observed the A-type star HD 162725, for which pre- 
vious measurements exist both with FORS 1 and later 
with ESPaDOnS at the Canada-France-Hawaii Telescope 
at much higher spectral resolution (_R=:65 000). Due to lack 
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Table 2. Targets discussed in this paper. 



Object 


V sin i 


log Mace 


IR 


Comp. 


log Toff 


log L 


M 


age 


name 


n — li 

[km s ] 


[Mo/yr] 


feature 




[K] 


[Lfi] 




[Myr] 


PDS2 


175" 






no 


3.98" 


1.68 


2.5 




HD 47839* 


120*" 






yes'='= 


4.57"" 


5.23"" 




r^O.3""" 


HD 95881 


50^" 


-8.04 


Si, PAH'" 


f f 

yggj J 


3.86"" 


1.84"" 


3.1 


>3'" 


HD 97048 


140'' 


-7.17 


PAH" 


no"® 


4.00" 


1.64" 


2.5 


>2'' 


HD 97300 




<-8.18 


PAH* 


h h 

yes"" 


4.02™™ 


1.54"'"' 


2.5 


>3"'"' 


HD 100453 


39'' 


-8.04 


PAH" 


yes" 


3.87" 


0.90" 


1.7 




HD 100546 


65^ 




Si, PAH 


no"® 


4.02" 


1.51" 


2.5 


>10'' 


HD 101412 


5^* 




Si, PAH 


yes"® 


3.98 


1.74^'" 


2.5 


rtf f f 


HD 135344B 


69*" 


-8.27 


Si^ PAH"" 


yes^-* 


3.82" 


0.91" 


1.6 


8±4''' 


HD 139614 


15' 


-7.99 


Si" 


yes'=" 


3.90" 


0.91" 


1.8 


>10'" 


HD 144432 


70' 


-7.07 


Si" 


yes'''' 


3.87" 


1.01" 


1.8 


10±5'" 


HD 144668 


100' 


-6.37 


Si'" 


yes" 


3.90" 


1.94" 


3.1 


0.6±0.4"'"' 


HD 150193 


100-'^ 


-7.29 


Si" 


yes'"'" 


3.95" 


1.38" 


2.2 


>2'' 


HD 152404 


18.5^ 




Si'' 


yes-* 


3.81" 


0.95" 


1.7 


~8"" 


HD 158643 


267'' /228'' 


-6.87 


Si" 


no"" 


4.00'* 


2.39'' 


3.9 


0.3"* 


HD 163296 


130' 


-7.12 


Si" 


no"® 


3.94" 


1.38" 


2.2 


5±2'" 


HD 169142 


66'' 


-7.40 


PAH"' 


yes"" 


3.91" 


1.16" 


2.0 


6±3'" 


VV Ser 


142° /229' 


-6.34 


Si" 


no"" 


3.95" 


1.27" 


2.0 


3.5±0.5""" 


HD 176386 


-8.11 


PAH"" 


yes'='= 


4.03 


1.69 


2.7 


>2'""' 


HD 179218 


60 


-6.59 


Si, PAH" 


no"" 


4.02" 


2.00" 


3.0 


1.3±0.5'" 


HD 190073 


12' 




Si, PAH"^ 


f f 
yes-* ■' 


3.95'' 


1.80'' 


2.7 


1.2±0.D 




St 


ars with debris disks 












HD 9672 


196" 




_bb 


no 


4.00"" 


1.42"" 


2.3 


~8"' 


HD 39060 


130'" 






no 


3.91"" 


0.95"" 


1.8 


r^l2-'-'-' 


HD 109573 


152'' 


<-8.53 


_bb 


yes*"*" 




1.32""" 


2.2 


8±2^^ 


HD 164249 








yes®® 


3.81 


0.48" 


1.4 


r^l2^-'-' 


HD 172555 








no 


3.89 


0.96" 


1.7 


~12^'^'^ 


HD 181327 


16" 






no 


3.81''"' 


0.54" 


1.4 


^12^^-' 


Ma 


gnetic Herbig 


Ae/Be star 


■s from previouE 


> studies 










HD 31648* 


go'' 




Si" 


no'''' 


3.94" 


1.14" 


2.0 




HD 104237* 


10'' 


-7.45 


Si" 


yes"" 


3.93" 


1.54" 


2.2 


2.0±0.5"'"' 


HD 200775* 


26*" 






yes" 


4.29" 


3.82" 




0.1±0.05"'"' 


V380 Ori 


200-'^ 




Si'" 


yes^^ 


3.97'* 




3.0 


1.5±0.5'" 


BFOri 


37' 




Si" 


no 


3.95" 


1.53" 


2.3 


r,in 




Standard Ap 


star with a 


weak magnetic 


field 










HD 162725 


32" 








3.99"'='= 


2.33""" 







Notes: The information on ?;sini values, mass-accretion rate, IR features, binarity, effective temperature, luminosity, and age was 
collected from the literature. The mass was deduced from the position of the star in the H-R diagram (see Sect. [5}. Objects 
marked with * are excluded from the correlation analysis (see Sects. El 14.31 and l5.ip . 

References: ""Vieira et al. ([2003|l . ''Bernacca & Perinotto (|1970|l . "Grady et al. (|1996|l . ''van den Ancker et al. lfT998ll . "Acke & 
■Waelkens p004)l . ^Hamidouche et al. p008|l . "this work, "Dunkin et al. (fT997)) . 'Hubrig et al. H2007b|) . ■' Alencar et al. ((2003]l . 
''Roye r et al . PPUT)) . 'Mora et al . PIKJT]) . "Slettebak (P^g^, " de la Reza & Pinzon p]OT), "Beskrovnaya et al. PIM)) , ^Ale cian 
et al. ((2008]) . " Lands treet et al. p008|l, ''B oersma et al. pOCT8|) , "Acke & van den A ncker (2004||, 'Siebenmorgenet al. (fT998)) . 
"Meeus et al. (gDSTJ, "Geers et al. (^007)) . "Sloan et al. (PDOg|) . '"van Boekel et al. (gDDSI,' "Przygodda et al. (HDDS}, 
"Pontoppidan et al. (|2007b|) . ""Siebenmorgen et al. ((2000)) . ''''Kessler-Silacci et al. (|2005|) . ""Okamoto et al. p004|l . ''''Schiitz et al. 
(|2005b|) . ""Dommanget & Nys (p94)) . •'^•'^Baines et al. ((2006)) . ""Corporon & Lagrange (fT999)) . ""Stelzer et al. (f2006)) . ''Chen et 

al. (|2006a|) . "Augcrcau et al. (^001]) . ''''Carmona et al. (PODT)) . " Steck lum et al. (HMSI, ""Reipurt h fc Z innecker dHMJ, 

""Roberge et a l. (,2002) , ""Grady et al. (pOOTll . ffStauffer et al. ((T995|) . ®®SIMBAD, ''''Eisner et al. ((2004)1 . ""Grady et al. ((2004|) . 
"Pirzkal et al. ((TWT)) . ""Schnerr et al. ((WT]) . ""Blondel et al. (d^, ""van den Ancker et al. ((TW)l . '"^Wade et al. (flOD7|, 
""Hughes et al. ((2008]) . ""di Folco et al. ((2004]) . """Debes et al. ((2008]) . ''"'Zuckerman & Song ((2004]) . """Folsom et al. pT)07]) . 
''''''Corcoran & Ray ((T998]) . """Walsh ((T980]) . •''•'^•'^ Wade et al. ((2005]) . """Pontoppidan et al. (|2007a|) . """Catala et al. ((2007]) . "'Thi 
et al. ((2IM]) . J^^Zuckerman et al. ((MIT]) . ''''''Simon et al. ((2000]) . and '"Natta et al. ((TOOT]) . 



of time we obtained only three series of two exposures, re- 
sulting in a magnetic field measurement (Bz) = 128±30G 
using the full spectrum and [B^) = 158±38G using hydro- 
gen lines. This star was previously observed by Bagnulo et 
al. (|2006p . who failed to detect a magnetic field in this star. 



The non-detection by these authors is probably due to the 
much lower spectral resolution in their spectropolarimetric 
observations with FOR S 1 (with a slit width of C.'S-l'.'O). 
Landstreet et al. ()2008p detected a weak magnetic field of 
the order of — lOOG with ESPaDOnS, and could also show 
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that the field is variable. We confirm their finding with 
our result, revealing a positive magnetic field of the order 
of 130-160 G, measured in metal and hydrogen lines (see 
Table[3]). We note, however, that an accurate ephemeris for 
HD 162725 is not known and presently it is not possible 
to prove whether our observations were carried out in the 
opposite rotational phase compared to the observations of 
Landstreet et al. (^005)) . 

Certainly, the advantage of using high resolution spcc- 
tropolarimeters such as ESPaDOnS at the Canada-France- 
Hawaii Telescope, and NARVAL at the Bernard Lyot 
Telescope at Pic du Midi Observatory (France) is indu- 
bitable, since higher spectral resolution observations pro- 
vide more detailed information about the behaviour of dif- 
ferent elements in the presence of a magnetic field and on 
the magnetic field topology. On the other hand, the use of 
hydrogen lines, which are usually not used for magnetic field 
measurements with ESPaDOnS and NARVAL, offers the 
unique opportunity to study the global structure of the de- 
tected magnetic field. The profiles of metal lines frequently 
exhibit conspicuous variations, which are signatures of the 
circumstcllar environment, stellar winds, of a non-uniform 
distribution of metals over the stellar surface, or of the 
presence of temperature spots (which are typical for late 
type stars). In our studies with FORSl, the longitudinal 
magnetic fields are measured in two ways: using only the 
absorption hydrogen Balmer lines or using the whole spec- 
trum including all available absorption lines, i.e. we use all 
lines, variable and non-variable, together. As hydrogen is 
expected to be homogeneously distributed over the stellar 
surface, the longitudinal magnetic field measurements sam- 
ple the magnetic field fairly uniformly over the observed 
hemisphere. 

3. Results 

The obtained magnetic field measurements are presented in 
Table [31 In the first two columns, we give the object name 
of the targets and the modified Julian dates of the middle 
of the exposures. In the third column we list the signal-to- 
noise ratio (SNR) calculated in the final one-dimensional 
spectrum around 4750 A. The measured mean longitudinal 
magnetic field {B^) using all absorption lines is presented in 
Col. 4. The measured mean longitudinal magnetic field {B^) 
using all hydrogen lines in absorption is listed in Col. 5. All 
quoted errors are la uncertainties. In Col. 6, we identify 
new detections by ND and confirmed detections by CD. 
We note that all claimed detections have a significance 
of at least Sct, determined from the formal uncertainties 
we derive. These measurements are indicated in bold face. 
In Figure [1] we show the significance distributions of our 
measurements. These distributions are different from the 
distribution shown by Wade et al. p007p for their sample 
of FORSl observations of Herbig Ae/Be stars. Assuming 
that only formal uncertainties apply, this suggests that the 
Herbig stars in our sample are probably weakly magnetic. 

Apart from the confirmed detections of a magnetic field 
in the stars HD 101412, H D 14466 8, and HD 190 073 (W ade 
et al. [20071 Hubrig et al. I2007b| Catala et al. [20071) . six 
other stars of our sample, PDS2, HD 97048, HD 100546, 
HD 135344, HD 150193, and HD 176386, show evidence for 
the presence of a weak magnetic field. About half of the 
stars with magnetic field detections possess longitudinal 
magnetic fields larger than 100 G. These stars are the best 





<B,> / 



Fig. 1. Distribution of the detection significance 
(i3z)/cr(B^) for our programme stars. Diagonal lines de- 
note the Herbig Ae/Be stars, while horizontal lines denote 
the debris disk stars. Top: Results using the full spectra 
(Col. 4 of Table [3]). Bottom: Results using only the hydro- 
gen Balmer lines (Col. 5 of Table [3]). 



candidates for future spectropolarimetric studies to analyze 
the behaviour of their magnetic fields over their rotational 
cycles to disclose the magnetic topology of their surfaces. 
For two Herbig Ae stars, HD 139614 and HD 144432, with 
previously detected weak magnetic fields at a significance 
level of 3cr (Hubrig et al. I2004bl Hubrig et al. I2007bp . the 
magnetic field in the present study was diagnosed only at a 
level of 1.8(7 and 1.6cr, respectively. Wade et al. (|2005p failed 
to detect a magnetic field in HD 139614 with quoted uncer- 
tainties of 25 G on two consecutive nights with the high res- 
olution spectropolarimeter ESPaDOnS. The marginal de- 



6 S. Hubrig et al.: Magnetic field measurements of Herbig Ae/Be stars with FORS 1/VLT 

Table 3. The mean longitudinal magnetic field measurements for our sample of Herbig Ae/Bc and debris disk stars 
observed with FORS 1 . 



Object 


MJD 


SNR 




— 7~D — ^ 


Comm. 


name 






[G] 


[G] 




PDS2 


54610.399 


1790 


103±29 


-4±70 


ND 


47839 


54609.968 


3375 


134±52 


141±62 




95881 


54609.993 


3745 


47±25 


24±20 




97048 


54609.137 


2980 


164±42 


188±47 


ND 


97300 


54609.047 


2495 


109±50 


96±53 




100453 


54610.022 


3820 


1±14 


-30±28 




100546 


54610.046 


4195 


89±26 


87±28 


ND 


101412 


54609.190 


3280 - 


312±32 


-454±42 


CD 




54610.081 


3525 - 


207±28 


-317±35 


CD 


135344B 


54609.243 


4020 


32±15 


41±30 






54610.145 


3970 


-38±11 


-37±12 


ND 


139614 


54610.201 


4010 


-40±25 


-58±32 




144432 


54609.383 


2595 


36±22 


32±42 




144668 


54610.238 


3710 


-62±18 


-92±27 


CD 


150193 


54609.092 


3300 


144±32 


— 252±48 


ND 


152404 


54609.275 


2815 - 


-107±37 


-53±23 




158643 


54609.275 


4165 


32±20 


14±23 




163296 


54610.255 


3535 


40±30 


105±40 




169142 


54610.175 


3875 


50±25 


23±32 




VVSer 


54610.332 


1050 


200±75 


376±141 




176386 


54610.272 


3230 - 


119±33 


121±35 


ND 


179218 


54609.360 


3855 


51±30 


28±34 




190073 


54609.410 


3900 


104±19 


120±30 


CD 


Stars with debris disks 


9672 


54345.389 


4365 


50±22 


42±32 




39060 


53455.193 


4205 


-16±14 


-36±22 




109573 


54610.438 


3815 


56±30 


41±32 




164249 


54610.301 


2970 


19±14 


1±25 




172555 


54610.287 


3875 


-35±20 


-31±28 




181327 


54610.364 


3495 


-38±16 


-75±26 






Standard Ap star with 


a weak maj 


jnetic field 




162725 


54549.413 


2435 


128±30 


158±38 





Notes: All quoted errors are Icr uncertainties. In Col. 6 we identify new detections by ND and confirmed detections by CD. We 
note that all claimed detections have a significance of at least 3(j, determined from the formal uncertainties we derive. These 
measurements are indicated in bold face. 



tections of magnetic fields in these stars during this ob- 
serving run can probably be explained by the strong depen- 
dence of the longitudinal magnetic field on the rotational 
aspect, i.e. on the rotation phase. For the Herbig Ae star 
HD 163296, we found no indication for the presence of a 
photospheric magnetic field, in agreement with our previ- 
ous studies (Hubrig et al. I2006bl Hubrig et al. l2007bp . 

No definite detection at a significance level of 3(t was 
achieved for any of the stars with debris disks. The only 
measurement close to the Scr level was obtained for the 
F5/F6V star HD 181327, which belongs to the (-12Myr 
old) /3Pictoris moving group, with a measured magnetic 
field (B^) = -75±26G. 

The star HD 101412, with the largest magnetic field 
strength measured in our sample stars, shows a change of 
the field strength by '^100 G during two consecutive nights. 
In Fig. [2] we present distinct Zeeman features detected at 
the positions of the hydrogen Balmer lines and the Can 
H&K lines. The H/3 line in the Stokes / spectrum is con- 
taminated by the presence of a variable emission in the line 
core and was not included in our measurements. 

Similar Zeeman features at the positions of the Can 
H&K lines were detected in four other Herbig Ae/Be 
stars, HD 139614, HD 144668, HD 152404, and HD 190073. 



In Fig. [3l we present the Stokes V spectra for these stars in 
the region around the Can doublet, together with our pre- 
vious observation of HD 190073, obtained with FORS 1 in 
2005. As we already reported in our earlier studies (Hubrig 
et al. I2004b[ [2006bl I2007b|) these lines are very hkely formed 
at the base of the stellar wind, as well as in the accretion 
gaseous flow, and frequently display multi-component com- 
plex structures in both the Stokes V and Stokes / spectra. 
In two Herbig Ae/Be stars, HD 31648 and HD 190073, such 
a structure was especially noticeable, and from their study 
we concluded that a magnetic field is present in both stars, 
but is most likely of circumstellar origin. For HD 31648, 
we detected a magnetic field (Sz) — 87±22 G. Using only 
the Ca II H&K lines for the measurement of circular polar- 
ization in HD 190073, we diagnosed a longitudinal magnetic 
field at a 2.8(t level, IB^ = 84±30 G. Our previous magnetic 
field measurement for HD 190073 of the order of 80 G was in 
full agreement with the high resolution spectropolarimet- 
ric measurements obtained for this star with ESPaDOnS 
(Catala et al. I2007p . who measured a longitudinal mag- 
netic field (Bz) = 74±10G using metal lines. The authors 
report that they could not see a Stokes V signal in the Call 
H&K lines, but the S/N ratio of their spectra is too low in 
that spectral region, resulting in a noise level of the order 
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Fig. 2. Stokes / and V spectra of the Herbig Ae/Be star 
HD 101412, with the largest detected magnetic field. Upper 
panel: Zeeman features in Hg, Hs, Can H&K, and He pro- 
files. Lower panel: Stokes / and V spectra in the vicinity of 
the H7 line. 

of 7x10^'^ in Stokes V per spectral bin of 0.0025 nm, which 
would not allow them to detect as weak a signal as seen in 
the LS D aver age, nor the level of signal reported by Hubrig 
et al. (|2006b|) . Our new observations reveal a rather large 
change (by ~0.5% of circular polarization) in the amplitude 
of the Zeeman features in the Can H&K lines compared to 
our observations from May 2005 (see the two upper spec- 
tra in Fig. [3]), indicating a conspicuous variability of the 
magnetic field. 

Wade et al. (|2007p mentioned the absence of Zeeman 
features in the Can H&K lines in their FORS 1 spectra 
of the Herbig stars HD 31648, HD 144432, and HD 144668. 
On the other hand, in the same study, they report that 
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Fig. 3. Stokes V spectra in the vicinity of the Can H&K 
lines of the Herbig Ae/Be stars HD 139614, HD 144668, 
HD 152404, and HD 190073. At the top we present our pre- 
vious observation of HD 190073, obtained in May 2005. The 
amplitude of the Zeeman features in the Ca 11 H&K lines ob- 
served in our recent measurement has decreased by ^^0.5% 
compared to the previous observations. 



their measurements were obtained using a variety of slit 
widths: 38 of 73 stars were measured with slit widths of 
0'.'8 to 1", and 35 of their 73 stars were measured with 
a slit width of 0'.'5. Using wide slits greatly degrades the 
resolution of the Stokes V spectra, possibly affecting the 
deduced longitudinal magnetic fields. Wade et al. (|2007p 
looked for this effect in their study, generally not finding 
strong evidence for it. 

4. Discussion of individual stars with magnetic field 
detections 

In the following we discuss the present knowledge of the 
various observed properties of stars with a magnetic field 
detection at a 3ct level, placing greater emphasis on the 
discussion of stars with stronger magnetic fields. Since a 
few stars of our sample have been previously observed by 
Wade et al. (|2007p . we also compare our results with the 
results obtained by these authors. 

4.1. Herbig Ae/Be stars with detected magnetic fields 

PDS 2: This star was observed with FORS 1 in November 
2004 by Wade et al. ((2007|) . A weak magnetic field of the 
order of 130-140 G was measured at a 3cr significance level 
using metal lines and the full spectrum, but could not be 
diagnosed from the measurements using hydrogen Balmer 
lines. We achieve a rather similar result: We detect a pos- 
itive magnetic field (Bz) = 103±29 G using the full spec- 
trum, but we are not able to detect the presence of a mag- 
netic field using hydrogen lines. For this star, no spatially 
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resolved observations of circunistellar matter or infrared 
spectral emission features were reported in the literature. 

HD 97048: No detection at a 3cr level was achieved 
by Wade et al. (|2007p . However, in our observations, 
the magnetic field in this star is rather strong. We 
measure (_Bz) ~ 161±42 G using the full spectrum and 
(Bz) = 188±47G using hydrogen lines. A large flaring disk 
with a mean disk inclination of 42.8° was reported by 
Doucet et al. (|2007p from observations with VISIR. No 
silicate emission band at 10 /im was found in the mid-IR 
spectra of HD 97048. The disk has been coronographically 
imaged by Doering et al. (|2007p . and its inclination angle 
appears to be closer to i=30° than the i=42.8° reported by 
Doucet et al. (I2007p . The disk has a radial surface bright- 
ness profile consistent with dust grain growth and settling, 
rather than a highly fiared geometry. No HH knots are vis- 
ible in the coronographic ACS F606W image presented by 
Doering et al. (|2007p . This is consistent with the at best 
low level of Mgii emission in the archival lUE data. No 
presence of a binary companion was detected in the study 
of Corporon & Lagrange (|1999p . The mid-IR spectrum is 
dominated by PAH emission, which arises mostly from an 
outer radius at 200-300 AU (van Boekel et al.'^DM). Both 
quiescent H2 emission in the near-IR (Bary et al. 2008J and 
rotational H2 emission of warm gas in the mid-IR (Martin- 
Zai'di et al. I2007P were reported. HD 97048 was detected 
in an XMM-Newton survey of Cha I. On basis of the nega- 
tive results of all searches for binarity, Stelzer et al. (|2004p 
ascribed the X-ray emission to the Herbig star. 

HD 100546: In this star we measure {B^} = 89±26 G us- 
ing the full spectrum and (B^) =87±28G using hydrogen 
lines. No detection at a 3cr level was achieved by Wade et al. 
(|2007p . Although no stellar companions have been reported, 
there seems to be evidence for a giant planet forming at 
~6.5AU (Acke & van den Ancker . The circumstel- 

lar disk was repeatedly resolved at optical, near-IR, mid-IR 
and mm wavelengths. Augereau et al. (|200ip derive a disk 
inclination of 51±3°, while the values from other authors 
are consistent within error ranges. There is a circumstel- 
lar envelope surrounding the disk (e.g. Grady et al. I200ip . 
Silicate emission and PAHs are seen in 10 /im spectra (e.g. 
Meeus et al. I200ip . Highly excited H2 gas was reported by 
Martin-Zai'di et al. (|2008p . who estimate that the warm H2 
is located within about 1.5 AU from the star. The star is de- 
tected in X-rays (Feigelson et al. [^0(151 Stelzer et al. 
with an Lx typical of Herbig Ae stars. No jet is seen in 
HST imagery in the optical (Grady et al. 120011 Ardila et 
al.HnnZD or in Lya (Grady et al. l2005ap . 

HD 101412: Wade et al. (PPHTP measured a positive 
magnetic field of the order of 500 G using hydrogen lines. 
During the first night we measured (B^) = — 312±32 G em- 
ploying the full spectrum and (i?z) = — 454±42 G on hy- 
drogen lines. The magnetic field was also detected on the 
second night: (Bz) = — 207±28G using the full spectrum 
and (Sz) = — 317±35 G on hydrogen lines. Interestingly, 
Wade et al. (|2007p found a negative magnetic field of the 
same order when they used metal lines for the measure- 
ments. Contrary to their results, our measurements show 
rather consistent results obtained using the full spectrum, 
metal lines and hydrogen lines, although the magnetic 
field from metal lines appears somewhat lower, probably 
due to CS contamination. From the measurements dur- 
ing the first night, using exclusively metal lines, we ob- 
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Fig. 4. SOFI Ks-band image of a potential faint companion 
to HD 101412. The spatial resolution is O'.'29/pix and the 
separation amounts to 3'.'3±0'.'2. 



tain (Bz) = -158±50G, and (B.,) = -95±46G for the sec- 
ond night. 

It is quite possible that HD 101412 belongs to a visual 
binary system. In Fig. [¥] we present a Ks-band image of 
HD 101412 obtained with SOFI at the NTT in La Silla. 
The data were taken in May 2000 in the framework of the 
ESO programme 65.I-0097(A) and were retrieved from the 
ESO science data archive. A faint candidate companion is 
detected at 3'.'3±0'.'2 separation to the southwest. The spa- 
tial resolution of the SOFI images was O'.'29/pix. Whether 
there is a real physical association between both stars re- 
mains to be confirmed. 

The very low wsini=7kms~^ suggests that we may 
see the disk nearly pole-on. Since this star exhibits the 
strongest magnetic field and hence is of special interest, we 
retrieved from the ESO archive three high-resolution UVES 
spectra of HD 101412 obtained in the framework of ESO 
programme 077.C-0521(A). An inspection of these spectra, 
recorded on three different dates, indicates variations in line 
intensities and line profiles. All lines seem to be variable, 
and we were not able to find any non-variable spectral line 
in the spectrum. The largest variations are observed be- 
tween the dates MJD 53872.540 and MJD 53920.494. A few 
examples of these variations are presented in Figs. [5] and [6l 
We tested a number of atmospheric models in the range: 
Teff =8000-11 000 K and log.g=4.0-4.5. The studied UVES 
spectra cover both Ha and H/3 Balmer lines. However, only 
the wings of the H/3 line can be used to obtain the at- 
mospheric parameters, as the Ha line exhibits strongly 
variable emission. We obtain as a best fit Toft=10000K, 
log 3=4.2-4.3 and a v sin i value of about 5 km s~^ . This re- 
sult is in good agreement with the study of Guimaraes et al. 
([20061) who published roff=10000±1000K, log.g=4.1±0.4, 
and 1; sini=7±l kms~^. A good fit is also obtained using 
the model reff=9000K and log 3=4.0, but the synthetic 
spectrum for these parameters shows many narrow lines 
which do not appear in the observed UVES spectra. The 
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Fig. 5. UVES spectra obtained on two different dates of 
the Herbig Ac/Bc star HD 101412 in the spectral regions 
around the Ha line (top) and around the H/3 line (bottom). 



spectrum is heavily contaminated by CS lines. Two such 
lines belonging to Fe i and Fe ii are clearly visible in Fig. [6l 
As we mentioned above, the low iisini value suggests 
that we observe HD 101412 close to pole-on, or that the 
star is rotating very slowly. The inclination angles of disks 
of Herbig Ae/Be stars (which are expected to be identical 
with the inclination angle of the stellar rotation axis) can 
be reliably derived only for resolved observations of disks. 
They are usually determined from millimeter observations, 
coronographic imagery, or near-IR interferometry. On the 
other hand, the orientation of the disk can be constrained 
using the emission profile shapes. Type I P Cygni profiles 
or single emissions all have i < 40°. Stars which alternate 
between type I and type III profiles have 40 < z < 55°, and 
stars with double emission all have i > 55°. This means 
that even for objects for which no good disk imagery ex- 
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Fig. 6. UVES spectra of the Herbig Ac/Bc star HD 101412 
in the spectral regions AA5526-5538A around the Cm, 
Mgi, and Fell lines (upper panel) and around the two Sill 
lines at A6347 (lower left panel) and A6371 (lower right 
panel). CS lines of Fei and Fell are clearly visible in the 
upper spectrum. Note the strong increase of the line inten- 
sities at MJD53920.494. 



ists, it is at least possible to constrain their inclination an- 
gles from the study of Mg ii spectral line profiles in the UV 
or from the line profile shape of Ha (e.g., Finkenzeller & 
Mundt \WMl Dunkin et al. [TM7|) . The shape of Ha in the 
UVES spectra of HD 101412 suggests an inchnation angle 
i > 55°, i.e. the star is viewed far from pole-on, in agree- 
ment with the recently published value 80° ±7° by Fedele et 
al. (PPDSI) who used VLTI/MIDI observations. The slow ro- 
tation could in principle be explained by binarity, where the 
components are synchronized with Pj-ot = ^orb- However, 
we detect neither spectral lines of the companion nor any 
variability of radial velocities in the UVES spectra. It can- 
not be excluded, though, that braking of the star's rotation 
is due to its rather strong magnetic field, in analogy with 
the slow rotation of magnetic Ap and Bp stars. Since the in- 
tensities of lines of different elements do not show opposite 
behaviour, which is usually observed in chemically pecu- 
liar Ap and Bp stars, we exclude the presence of chemical 
spots on the stellar surface. Also the study of temperature 
sensitive and insensitive lines indicates the absence of tem- 
perature spots. Therefore we conclude that the observed 
spectrum variability has a CS origin. Silicate emission and 
lower amounts of PAHs are seen in mid-IR spectra (e.g., 
van Boekel et aL l^DUSl Geers et al. lMTT)) . w hile th e SED is 
consistent with a flat disk type (Geers et al. I2007p . 

HD 1353UB=SAO 206462: The detected magnetic field 
in this star is very weak. We observe a change of the polar- 
ity between the first and the second observing night, but 
the measurement on the first night, using the full spec- 
trum, resulted in a strength of the longitudinal magnetic 
field (i?z)=32±15G which is only at a 2.1(7 significance 
level. On the second night, the magnetic field was detected 
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at a 3.2a level using the full spectrum. (Bz) = — 38±12G. 
As is mentioned in the literature, there has been frequent 
confusion between HD 135344 and the ^20" close, nearby 
IR so urce SAO 206462, a mid-F star (Coulson & Wahher 
[USB, also referred to as HD 135344B. Thus SAO 206462 
was frequently improperly called HD 135344. A further bi- 
nary pair was found at '^5'.'8 separation from SAO 206462 
by Augereau et al. (|200ip . Mid-IR spectra show silicate ab- 
sorption around 10 ^m (e.g. Acke & van den Ancker [2004p 
and PAHs (e.g. Sloan et al. [MIS]) . Dent et al. (^005)) de- 
rived a disk inclination of 11°, obtained from measurements 
of CO gas. The presence of warm H2 gas is reported by Thi 
et al. (|200ip . which should be treated with caution due to 
the low S/N and low spatial resolution of ISO data. 

This star has recently been identified as the host of 
a transitional or pre-transitional disk (Brown et al. I2007| 
Pontoppidan et al. l2008p . The disk is coronographically de- 
tected at 1.1 fiui and 1.6/zm (Grady et al., in preparation) 
with a radial surface brightness profile consistent with grain 
growth and settling. No jet is seen in HST optical corono- 
graphic imagery (Grady et al. I2005bp . FUV excess light 
and emission were detected by FUSE, with a mass-accretion 
rate reconcilable with the Br7 measurements (Garcia Lopez 
[20061 Grady et al, in preparation). The bulk of the FUV 
emission is consistent with stellar activity, rather than with 
accretion. 

HD 144668=HR 5999: In our previous study (Hubrig et 
al. I2007ap . we showed that the magnetic field of this star 
is variable, with the strength varying between —75 G and 
-1-166 G. Strong Zeeman features are clearly visible at the 
position of the Can H&K lines and hydrogen Balmer lines, 
similar to those detected in our previous observations (see 
Fig. [3|). We confirm the presence of the magnetic field, 
measuring = — 62±18 G using the full spectrum and 
(i3z) = — 92±27G using hydrogen lines. In stellar applica- 
tions, the Stokes V profiles present disk-integrated observa- 
tions. Generally, stellar magnetic fields are not symmetric 
relative to the rotation axis, so that the polarization signal 
changes with the same period as the stellar rotation. The 
most simple modeling includes a magnetic field approxi- 
mated by a dipole with an axis inclined to the rotation axis. 
Recently, Auriere et al. (|2007p presented measurements of 
weak stellar magnetic fields using the cross-correlation tech- 
nique "Least Squares Deconvolution" , originally introduced 
by Donati et al. (|1997p . They showed that the presence of 
a weak magnetic field can be diagnosed from the detection 
of Stokes V Zeeman signatures, even if the longitudinal 
magnetic field is very small, as is usually observed dur- 
ing a cross-over phase. In our observations of a few Herbig 
Ae/Be stars, even at the very low resolution of FORSl, 
similar Stokes V Zeeman features are observable in strong 
Can H&K and hydrogen Balmer lines, indicating the pres- 
ence of a cross-over phase. In Fig. [7] we show an example of 
Stokes V line profiles of H7 and H/3 in HD 144668, confirm- 
ing the presence of a magnetic field in this star at a rotation 
phase where the longitudinal magnetic field is rather small. 

HD 144668 belongs to a quadruple system (Dommanget 
& Nvs ll994)) . although the orbits and physical relations be- 
tween the companions are not certain. The nearest neigh- 
bour at a separation of 1'.'4 appears to be a TTauri star 
(Stecklum et al. I1995p . A broad but rather shallow 10 /im 
silicate emission feature indicates grain growth (van Boekel 
et al. l2005p . Preibisch et al. Ij2006p derived an inclination of 
58° from MIDI/VLTI observations and found that the disk 
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Fig. 7. Stokes / and V spectra of the Herbig Ae/Be star 
HD 144668. Upper panel: Zeeman feature in the H7 line. 
Lower panel: Zeeman feature in the H/3 line. 



is very compact with an outer edge of '^2.6 AU, potentially 
caused by a further, even closer, yet undetected companion. 
Recently, HD 144668 was for the first time resolved from its 
TTauri companion in X-rays (Stelzer et al. I2009p . This ob- 
servation showed that the bulk of the X-ray emission comes 
from the companion, but a weak source is associated with 
HD 144668 (or the hypothesized sub-arcsecond companion). 

HD 150193: No attempt had been made to measure the 
magnetic field in this object. The detected magnetic field 
is comparatively strong: we measure (_Bz) = — 144±32G us- 
ing the full spectrum and {B^) = — 252±48G using hydro- 
gen lines. A binary companion at I'.'l distance (Reipurth 
& Zinnecker ,19931) is a T Tauri star (e.g. Carmona et al. 
I2007p . Fukagawa et al. (j2003p resolved the circumstellar 
disk in the near-IR and derived a disk inclination of 38±9°. 
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Broad 10 fim silicate emission is seen in mid-IR spectra (e.g. 
Ackc & van den Ancker I2004p . Warm H2 gas is reported 
by Thi ct al. ||200T|) . Similar to the case of HD 144668, re- 
solved Chandra images show that most of the X-ray emis- 
sion comes from the T Tauri companion, but weak emission 
is attributed to the Herbig star (Stelzer et al. I2006p . 

HD 176386: A magnetic field is definitely present: we 
measure (Sz) = — 119±33G using the full spectrum and 
(i?z) = — 121±35 G using hydrogen lines. This star is part 
of a triple system (Dommanget & Nvs 11991)) . Grady et al. 
(|1993p reported evidence for accretion of circumstellar gas 
onto the star. Martin-Zaidi et al. (|2008p detected warm H2 
gas. This source shows the highest H2 column density in 
their entire sample. Both IR (Siebenmorgen et al. l2000p and 
far-UV spectroscopy (Martin-Zaidi et al. I2008P suggest that 
the circumstellar matter has the form of an envelope rather 
than the form of a disk. PAHs are detected in mid-IR spec- 
tra (e.g. Siebenmorgen et al. l2000p . Surprisingly, HD 176386 
is one of the few Herbig Ae/Be stars confirmed to be X-ray 
dark (log Lx[erg/s] < 28.6). Stelzer et al. pM)) showed on 
basis of Chandra imaging that apparent X-ray detections 
of this star in previous lower-resolution observations are to 
be attributed to the companion at 4" separation. 

HD 190073: A magnetic field of (Sz)==74±10G was 
detected using photospheric metal lines with ESPaDOnS 
(Catala et al. I2007p . Our old measurements (Hubrig et al. 
I2006ap revealed a magnetic field of the order of 80 G at 
the 2.8(7 level, in full agreement with the high resolution 
ESPaDOnS spectropolarimetric data, whereas Wade et al. 
(|2007p failed to detect the magnetic field in both previ- 
ous measurements carried out with FORS 1 in November 
2004. Our new measurements = 104±19 G using the 
full spectrum and (_Bz) = 120±32 G using hydrogen lines 
confirm the results presented by Catala ct al. (|2007p . Baines 
et al. (|2006p reported the presence of a potential com- 
panion. Silicates and PAHs are seen in the mid-IR (e.g. 
Boersma et al. I2008p . Eisner et al. (|2004p showed that in- 
terferometric visibilities are consistent with a disk close to 
face-on, although a significant nonzero inclination cannot 
be ruled out. 

4.2. Stars with measurements suggestive of tlie presence of 
magnetic fields 

HD 47839=15 Mon: This is a massive O star classified as 
a pre- main sequence star in SIMBAD. According to Kaper 
(|1996p and Walborn (|2006p it shows distinct peculiarities 
in the spectra, which could be typical for stars possessing 
magnetic fields. The measured magnetic field is of positive 
polarity at a 2.6a significance level, using the full spectrum. 
The star should be considered as an important target for 
follow-up studies of magnetic fields in massive stars. 

HD 97300: Even though the measured field is only at 
the ~2(T level, we detect a few distinct Zeeman features in 
the Stokes V spectrum, supporting the evidence of a mag- 
netic field. This star is a good candidate for future magnetic 
field measurements. No quiescent H2 gas was found (Bary 
et al. I2008p . Siebenmorgen et al. (|1998p reported the pres- 
ence of a large circumstellar elliptical ring of PAH emission 
with thousands of AU in diameter and an indication for 
the possible existence of a candidate companion. It was 
also detected as an X-ray source in ROSAT and Chandra 
observations (Hamaguchi et al. 120051 Stelzer et al. I2006P 



HD 139614: Our previous measurements of HD 139614 
revealed the presence of a weak magnetic field in the 
range from -116 G to -450 G. (Hubrig ct al. l2004b[[2006bl 
I2007bp . The current diagnosis of magnetic fields in this star 
yielded no detection at a 3tT level, neither in Balmer nor 
in metal lines, with error bars as low as 25 G. In analogy 
with our previous detections, the measured magnetic field 
is of negative polarity and Can H&K lines in HD 139614 
show distinct Zeeman features. The nearly featureless mid- 
IR continuum of HD 139614 shows small amounts of sili- 
cates, but no PAH (Keller et al. lMIH)) . 

HD 144432: For this object we measured in the past a 
magnetic field in the range from — 94G to — 119G (Hubrig 
et al. l2004bl[2006bll2007bl) . Dent ct al. mm found CO gas 
and derived a nearly face-on disk orientation. HD 144432 is 
a binary with ~1'.'4 separation (Dommanget & Nvs ll994|) . 
Carmona et al. (|2007p reported the detection of Ha emis- 
sion and Li absorption in the K5Ve T Tauri companion. 
The companion is an X-ray source, and weaker X-ray emis- 
sion is also detected from the position of the Herbig Ae/Be 
star (Stelzer et al. I2009p . The 10 /im silicate emission band 
shows rather pristine dust, while the SED points towards a 
flat disk, as reported by Meeus et al. (|200ip . 

HD 152404=AK Sco: The magnetic field 
(Sz) = 107±37G is detected at the 2.9cr level using 
the full spectrum. The Stokes V spectrum exhibits distinct 
Zeeman features at the position of higher number Balmer 
lines and of the Can H&K lines (see for example Fig. 
This star is one of the best candidates for future spec- 
tropolarimetric observations. The Chandra image presents 
one faint X-ray source associated with HD 152404 (Stelzer 
et al. I2006P . Some authors refer to this target as a T Tauri 
star, as its spectral type may be at the borderline. The 
10 /xm silicate emission feature indicates very pristine dust 
(e.g. Przygodda ct al. [2003)) . 

VVSer: The magnetic field is diagnosed only at a sig- 
nificance level of 2.7(7, but additional measurements are de- 
sirable. Alonso-Albi et al. (|2008p carried out observations 
at millimeter and centimeter wavelengths towards VV Ser 
using the Plateau dc Burc Interferometer and the Very 
Large Array to compute the SED from the near infrared 
to centimeter wavelengths. The modeling of the full SED 
has provided insights into the dust properties and a more 
accurate value for the disk mass. The mass of dust in the 
disk around VVSer was found to be about 4x10~^Mq, 
i.e. 400 times larger than previous estimates. Due to the 
faintness of the star (my = 11. 6), only two series of two sub- 
exposures were taken, each one with an exposure time of 
15 min. Pontoppidan et al. (j2007bp compared Spitzer in- 
frared data with a disk model and found that the disk is 
nearly edge-on and self-shadowed by a puffed-up inner disk 
rim. The 10 /zm silicate emission feature is rather weak. 

4.3. The observed circumstellar properties of Herbig Ae/Be 
stars with magnetic field detections discussed in previous 
studies 

HD 3I648: Based on interferometric observations, Eisner et 
al. (|2004p determined a disk inclination of ~30° and re- 
ported that a binary model can be ruled out with a high de- 
gree of confidence. Further, Eisner (|2007p found a visibility 
increase in the Br7 line for this star, as expected for mag- 
netospheric accretion. Pietu et al. (|2007p had obtained an 
inclination angle of ~35° from mm measurements. Silicate 
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emission was detected in the mid-IR around 10 /im (e.g. 
Ackc & van den Ancker I2004p , while the presence of warm 
H2 gas is reported by Thi et al. (|200ip . 

HD 104237: This star is at least a quintuple system 
(Grady et al. I2004p . where some components are TTauri 
stars. A further spectroscopic companion is suggested by 
observations of Bohm et al. (|2004p and Baines et al. (|2006p . 
The Herbig star itself is the brightest X-ray source of the 
group (Stelz er et al. The disk is inchned by -18° 

(Grady et al. l2004p . Silicate emission is seen in mid-IR spec- 
tra (e.g. Acke & van den Ancker [2004p . Martin-Zaidi et al. 
(|2008p found excited and hot H2 gas from far-UV spec- 
tra, while TatuUi et al. (|2007p explain the origin of the Br7 
emission by an inner disk wind, originating at about 0.5 AU 
from the star. 

HD 200775: This object has a companion detected at 
2'.'25 (Pirzkal et al. I1997p . while spectroscopic observa- 
tions suggested further binarity of the main component 
(e.g. Miroshnichenko et al. ll998p . which was repeatedly dis- 
cussed in the literature and later confirmed by Monnier et 
al. (|2006p using interferometry. 

V380 Ori: This source is actually a triple system with 
separations of S'.'O and S'.'O (Dommanget & Nvs ll994)) . while 
the primary itself has a companion at 0'.'15 (Leinert et al. 
I1997P which is probably another Herbig Ae/Be star. Even 
the superior spatial resolution of Chandra is not able to 
resolve this binary, such that the origin of the X-ray source 
remains obscure (Stelzer et al. l2006p . The broad but shallow 
10 fim silicate feature indicates grain growth, as confirmed 
by van Boekel et al. ^UU^. 

BF Ori: As Mora et al. (|2004p showed, the spectra ex- 
hibit circumstellar line absorptions with remarkable varia- 
tions in their strength and dynamical properties. High ve- 
locity gas is observed simultaneously in the Balmer and 
metallic lines and the gaseous circumstellar environment 
of this star is very complex and active. Silicate emission 
is detected in 10 /im spectra (e.g. Acke & van den Ancker 

mm . 

5. Searching for a link between the presence of a 
magnetic field and fundamental stellar 
characteristics 

Vink et al. (|2002p presented their results of Ha spectropo- 
larimetric observations of a sample of 23 Herbig Ae/Be 
stars, pointing out the possibility of the existence of a phys- 
ical transition region in the H-R diagram from magneto- 
spheric accretion, similar to that of cTTS, at spectral type 
A to disk accretion at spectral type B. The main differ- 
ence between these scenarios is that in the former case the 
stellar magnetic field truncates the accretion disk at a few 
stellar radii and gas accretes along magnetic channels from 
the protoplanetary disk to the star, while in the latter case 
the accretion flow is not disrupted by the field. Our sample 
consists of 21 Herbig Ae/Be stars of spectral classification 
B9 and later spectral types and six debris disk stars. Since 
the observations of the disk properties of intermediate mass 
Herbig stars suggest a close parallel to cTTS, it is quite pos- 
sible that magnetic fields play a crucial role in controlling 
accretion onto and winds from Herbig Ae stars, similar to 
the case of the lower- mass cTTS. Evidence for disk accre- 
tion in Herbig stars from optical emission line profiles was 
presented by MuzeroUe et al. (|2004p . However, contrary to 




LogTrf, [K] 

Fig. 8. H-R diagram for all Herbig Ae (filled circles) and de- 
bris disk (filled squares) stars from Table[2]on the pre-main 
sequence models from Siess et al. (|2000p . Dotted curves 
show isomass evolutionary tracks and solid jagged contours 
indicate the size of the convective envelopes labeled in units 
of the stellar radius. 



the advances achieved in magnetic studies of cTTS, there 
is still no substantial observational evidence demonstrat- 
ing the strength, extent, and geometry of magnetic fields 
in Herbig Ae stars. We are aware of the fact that our ob- 
servations do not present a systematic monitoring of the 
magnetic fields of Herbig stars over the rotation period and 
are just snapshot observations of (Bz) values over two vis- 
itor nights. However, the magnetic field measurements in 
these stars are rare due to the very small number of spec- 
tropolarimetric facilities on large telescopes, and presently 
no other magnetic field data are available. Still, as we show 
in the next sub-sections, a few hints and trends can be es- 
tablished with the obtained data. The search for a link with 
other stellar properties is important to put preliminary con- 
straints on the mechanism responsible for magnetospheric 
activity. As we show in the next sub-sections, we establish 
for the first time preliminary trends. 

The H-R diagram for all Herbig Ae and debris disk stars 
from Table [2] is shown in Fig. [H In the following, we study 
the disk properties, binarity, age, and X-ray emission for 
the Herbig Ae stars and debris disk stars. We include in the 
definition of 'Herbig Ae' stars targets with spectral type B9 
to mid-F. We have estimated the individual stellar masses 
from an interpolation of the evolutionary tracks from Siess 
et al. (pOgPl) . and list the derived values in Col. 8 of Tabled 
As can be seen from Fig. |S1 the bulk of our sample has 
masses between 2 — 3 Mq . A comparison with evolutionary 
calculations shows that the majority of the stars in our 
sample can be considered fully radiative: the depth of the 
convection zone is less than 1 % (cf. contours in Fig.[5|). 

In the remainder of this section, we discuss our search 
for dependencies between the magnetic field and other char- 
acteristics of our targets. We have examined various pairs 
of relevant parameters with correlation tests implemented 
in ASURV (Astronomy Survival Analysis Package; Lavalley 
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Table 4. Results from statistical analysis of parameter 
pairs. 



Table 5. Longitudinal magnetic fields of four previously 
studied Herbig Ac stars with detections at the 3a signifi- 
cance level. 



Param 1 


Param 2 


iV 


T"Kcnd 


pSpcar 


log Mace 


Age 


16 


0.05 


0.08 




log Mace 


16 


0.68 


0.57 


log Lx 


Age 


19 


0.13 


0.28 




log Lx 


19 


0.09 


0.07 


Bz 


log (Lx/ibol) 


19 


0.14 


0.19 


Bz 


^rot 


20 


0.50 


0.42 




Age 


27 


0.03 


0.05 



Notes: Parameter pairs are given in Cols. 1 and 2. A'^ is the 
number of data points, TKond and pspoar denote the probability 
for no correlation according to Kendall's and Spearman's rank 
order correlation test, i.e. small numbers indicate that a 
correlation is present. 
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Fig. 9. Mass-accretion rate versus age. Open circles denote 
Herbig Ae stars and open squares indicate debris disk stars. 



et al. ll992p . The probabilities for a correlation between two 
given parameters are summarized in Tabled 

5.1. Accretion rate 

The recent results of low resolution linear spectropolari- 
metric observations of Herbig stars in Ha, H/3, and H7 
by Mottram et al. (|2007p support the presence of mag- 
netospheric accretion in Herbig Ae stars. The authors de- 
tected intrinsic line-polarization signatures suggesting that 
the magnetic accretion scenario generally considered for 
cTTS may be extended to Herbig Ae stars, but that it 
may not be extended to early Herbig Be stars, for which 
the available data are consistent with disk accretion. 

The number of Herbig Ae stars with measured mass- 
accretion rates is still rather low (see Col. 3 in Table [5]). 
In Fig. [21 we show the dependence of mass-accretion rate 
on age for the stars in our sample. Obviously, the highest 
mass- accretion rates tend to be found in younger stars. 



Object 




name 


[G] 


HD 139614 


-93±14 


HD 144432 


-111±16 


V380 Ori 


460±70 


BFOri 


-144±21 



In the following figures we use for each studied Herbig 
star the maximum measured value of the longitudinal field 
measured from hydrogen lines without taking into account 
the polarity of the field. In addition, we increased the sam- 
ple of magnetic Herbig Ae stars by adding to our sample 
four stars with previous detections of photospheric mag- 
netic fields, namely HD 139614, HD 144432, V380Ori, and 
BF Ori. As mentioned in Sect.[Tl our previous observations 
of HD 139614 and HD 144432 revealed the presence of a 
photospheric magnetic field of the order of ~100G (Hubrig 
et al. I2007bl) . The longitudinal magnetic field strengths of 
HD 139614, HD 144432, V380Ori, and BFOri are listed in 
Table El 

The values of Mace were derived by Garcia Lopez et al. 
(|2006p from the measured luminosity of the Br7 emission 
line, using the correlation between L(Br7) and the accre- 
tion luminosity Lace, established by MuzeroUe et al. (|1998p 
and Calvet et al. (|2004p . The correlation used is empiri- 
cal, which makes no assumptions on the origin of Br7. An 
important result found recently by Kraus et al. (|2008p is 
that the Br7 line can trace both mass infall and outflow, 
implying that Br7 is probably only an indirect tracer of the 
mass-accretion rate. The authors used the VLTI/ AMBER 
instrument to spatially and spectrally resolve the inner (<5 
AU) environment of flve Herbig Ae/Be stars (HD 163296, 
HD 104237, HD 98922, MWC 297, and V921 Sco) in the Br7 
emission line as well as in the adjacent continuum. The 
quantitative analysis for HD 98922 reveals that the line- 
emitting region is compact enough to be consistent with 
the magnetospheric accretion scenario, and for the stars 
HD 163296, HD 104237, MWC 297, and V921 Sco the au- 
thors identify an extended stellar wind or a disk wind as 
the most likely line-emitting mechanism. We have not yet 
searched for a magnetic field in the star HD 98922, but we 
observed the star HD 163296 twice. Both observations were 
non-detections, supporting the results of Kraus et al. (|2008p 
that the magnetospheric accretion scenario does not work 
for this star. No significant magnetic field detection was 
achieved for HD 104237 by other authors (e.g. Wade et al. 
I2007p . In Fig.[Tni we show the correlation between mass ac- 
cretion rate and measured longitudinal magnetic fields for 
our Herbig Ae sample. While we do not see a simple cor- 
relation between the magnetic field strength and the mass 
accretion rate in the data, the observed values are in the 
range of predictions from magnetospheric accretion models. 

Magnetospheric accretion models describe the interac- 
tion between a dipolar stellar magnetic field and a sur- 
rounding accretion disk assuming pressure equilibrium. The 
analytical approach yields equations that relate the mag- 
netic field strength to the system parameters (mass, radius, 
accretion rate, and rotation period; see e.g., Koenigl 119911 
Shu et al. ^MM- Johns-KruU plilT)) has shown that for 
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the case of cTTS, the models developed by different inves- 
tigators make consistent predictions on the magnetic field 
strength. 

We examine the magnetospheric accretion scenario for 
our sample of Herbig Ae stars using the expression for 
the surface equilibrium field given by Koenigl (|199ip . We 
present this relation in the form adopted by Johns-KruU et 
al. (figMl) : 



5/6 



d 



(5) 



where -Beq denotes the strength of the equatorial magnetic 
field in kG, the parameter e expresses the stellar rotation 
rate in terms of rotation at the inner boundary of a trun- 
cated accretion disk, the parameter (3 describes how effec- 
tively the stellar magnetic field couples to the inner regions 
of the disk, denotes the stellar mass in solar masses, 
M the accretion rate in 10~^ solar masses per year, i?, 
the stellar radius in solar radii, and the stellar rotation 
period in days. 

In Fig. [ini the dotted line represents the prediction from 
Eq. [5]for a canonical star with = 2.5 Mq, i?* = 2.5 R©, 
and ~ 0.5 d. We find that the observed field strengths 
qualitatively support the magnetospheric accretion model 
for Herbig Ae stars, although the correlation analysis does 
not yield a positive result for (B^) vs. Mace (Tabled]). 

It must be kept in mind that Eq. [5] refers to the equa- 
torial dipolc field, while our measured values represent an 
average of the longitudinal component over the stellar sur- 
face and they also depend on the viewing angle. 

Future studies of the magnetic field topology and an 
improvement of indicators of mass-accretion rates are im- 
portant to understand the role of magnetic fields in the 
dynamics of the accretion processes in Herbig Ae stars. 

5.2. X-ray emission 

In Fig.[8]it is demonstrated that our targets have very shal- 
low or completely absent convection zones. Standard dy- 
namo theory does not predict magnetic field generation for 
fully radiative stars. Our targets are also not hot enough to 
drive strong radiative winds. Therefore, any X-ray activity 
is expected to decay with the dissipation of the primordial 
field. Nevertheless, many Herbig Ae/Be stars are known to 
be X-ray sources (Zinnecker & Preibisch 119941 Hamaguchi 
et al lMISl Stelzer et al. lMISllMjgi and references therein). 
In the absence of a theoretical model for X-ray production 
in intermediate-mass stars, the detections have often been 
ascribed to known or assumed late-type T Tauri star com- 
panions. However, recent Chandra imaging studies have re- 
solved Herbig Ae/Be stars from most known companions 
and still came up with very high detection rates (Stelzer et 
al. 120061 120091) . 

We collected X-ray luminosities from the literature, 
adopting for each star the value from the instrument with 
the highest spatial resolution available. The logLx and 
log Lboi values and the corresponding references are pre- 
sented in Table [6| Many of the stars of our sample were 
included in dedicated X-ray imaging studies with Chandra 
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Fig. 10. The strength of the longitudinal magnetic field as 
a function of accretion rate. Overplotted is the relation pre- 
dicted by magnetospheric accretion models (see Eq. [5]) for 
assumed stellar parameters of Af, = 2.5 M©, i?* = 2.5 R©, 
and P* = 0.5 d. Filled circles denote Herbig Ae stars with 
a 3f7 magnetic field detection, while open circles denote 
Herbig Ae stars with non-detections. Squares denote stars 
with debris disks, of which none has a Scr magnetic field 
detection. 



aimed at resolving them from their visual late-type com- 
panion stars (see references in Table [S]). Consequently, only 
a few of our targets have known unresolved sub-arcsecond 
or spectroscopic companions that might be responsible 
for the observed X-ray emission. In particular, AK Sco 
(=HD 152404) is a spectroscopic binary (Andersen et al. 
I1989p . V380 Ori is both a visual (0'.'154) and a spectroscopic 
binary (Leinert et al. 119971 Corporon & Lagrange I1999p . 
and HR479 6 (=HD 109573) is a 7'.'7 binary (Jayawardhana 
et al.HnnSD that could not be resolved with ROSAT. 

In Fig. [Tl] we plot logLx over the age for our sample 
stars. It is impossible to exclude categorically that there 
are further, as yet unknown, companion stars responsible 
for the X-ray emission. However, given the high detection 
rate and high spatial resolution of the X-ray studies, this is 
unlikely to be the case for all stars in this sample. 

Tout & Pringle p995p have suggested a mechanism that 
can give rise to intrinsic X-ray emission from Herbig Ae/Be 
stars. In their model, a dynamo can be sustained in a ra- 
diative star by rotational shear. As a result of the decrease 
of the available rotational energy with time, the X-ray lu- 
minosity behaves as 



Lx — Lx.Q • I 1 + — 
To 



(6) 



i.e. it decays from an initial value of Px.o given by the 
shear energy. Tout & Pringle (|1995p showed that an X-ray 
luminosity of about 1% of the bolometric luminosity can be 
maintained for a duration of '^l Myr, and thereafter decays 
rapidly according to the law presented in Eq. [Sj The curves 
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Table 6. X-ray emission observed in Herbig Ae stars and 
debris disk stars. 



Object 


log Lx 


Instrument 


1 / r It \ 

log (Lx/I/boi) 


name 


[erg/ s] 






HD 97048 


29.58 


XMht 


-5.64 


HD 97300 


29.31 


Chandra' 


-5.81 


tiD i(JU45o 


28.82 


Lhanara 


—5.66 


tiu i(JU54D 


28.93 


Lhanara 


—6.16 


tlU io0o44J3 




KUoAl 


^ A SO 




zo.y 


Chandra^ 


RQ 




28.3 


V_>/ bUil bUil Lb 


—7.22 




29.32 


\^ 1 bai bai Lb 


— ^ fid 


HD 152404 


29.09 


Chandra^ 


—5.44 








<-_f; qq 




9Q R 


^1 bat bat a 


— o.ou 


HD 169142 


29.08 


Chandra*^ 


-5.66 


HD 176386 


<28.56 


Chandra!' 


<-6.71 


HD 39060 


26.5 


Chandra'' 


-8.03 


HR 109573 


29.4 


ROSAT> 


-5.50 


HD 172555 


<26.9 


Chandra'' 


<-7.64 


HD 181327 


<29.4 


ROSAT 


<-4.72 


V380 Ori 


30.96 


Chandra'' 


-4.66 


BF Ori 


<30.2 


ROSA'f 


<-4.91 



Notes: The log Lx values have been collected from the 
literature, while the log Lboi are from Table [2] 
References: "Stelzer et al. (f2004ll . ''Stelzer et al. (f2006)) . 
'^Collins et al. (|2009|l . ''Grady et al. ApJ, submitted, ''Stelzer et 
al. (f2009)l . •''Berghoefer et al. (fT996)l . ^Swartz et al. p005|) . 
''Grady et al. (|2007p . 'Robrade, priv. comm. ■'Stelzer & 



Ne uhause r (pOOOjl . 
al. (|2005)) . 



'^Feigelson et al. (|2006|l . and 'Hamaguchi et 



plotted in Fig. [TT] represent the predicted decay of the X- 
ray emission for to = 1 Myr and 2 Myr. The initial X-ray 
luminosity Lx,o was estimated from Eqn. 4.4 of Tout & 
Pringle (IT995]) for an assumed = 2.5 M©, i?* = 2.5 M©, 
and Lboi = 50 L© (dash-dotted line). Two other solutions 
with a factor of two lower Lx.o £^rG also shown. 

No clear dependence of the X-ray emission level on the 
age is noticeable in the data but the absence of stars with 
high X-ray luminosities at an advanced age is in agreement 
with the expected decay of the shear dynamo. 

In Figs . [T^ and [T51 we present the strength of the mag- 
netic field plotted versus logLx and log(Lx/iboi)- In both 
figures we find a hint for an increase of the magnetic field 
strength with the level of the X-ray emission, also supported 
by the correlation analysis which yields >90% probabil- 
ity that a correlation is present between {B^} and logLx 
(Table 3]). This could suggest a dynamo mechanism respon- 
sible for the coronal activity in Herbig Ae stars. On the 
other hand, we should keep in mind that the star with the 
strongest magnetic field in both figures is the spectroscopic 
binary V380Ori with (B^) = 460±70G, and it is not clear 
yet whether the X-ray emission originates from the primary 
or from the companion. The filled circle with the lowest 
strength of the longitudinal magnetic field, (B^) = 37±12 G, 
belongs in both figures to the Herbig Ae star HD 135344B. 

Pcvtsov et al. (|2003p have derived a universal relation 
between magnetic flux and X-ray luminosity that holds over 
10 orders of magnitudes from the quiet Sun, over solar ac- 
tive regions to active late- type field dwarfs. The unique 




6 8 10 
Age [Myr] 



Fig. 11. X-ray luminosity versus age. Symbols are identical 
to Fig. [HI The decay of Lx predicted by the shear dynamo 
of Tout & Pringle (|1995p is shown for different values of 
the initial parameters to and Lx,o- The dotted line repre- 
sents io=lMyr and Lx,o=5 x 10"^"^, the dashed hne repre- 
sents to=2 Myr and Lx,o=5 x 10^^, and the dash-dotted line 
io=lMyr and Lx,o=9 x lO^^. 



power-law relation Lx ~ ^f^-^^ was taken as evidence that 
the coronae of the Sun and the stars arc heated by the 
same kind of structures. We test whether Herbig Ae stars 
obey the same relation. The magnetic fluxes of our stars 
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Fig. 12. The strength of the longitudinal magnetic field 
plotted over the X-ray luminosity. Symbols arc identical to 
Fig.IlOl 
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Fig. 13. The strength of the longitudinal magnetic field 
plotted over the ratio between X-ray luminosity and bolo- 
metric luminosity. Symbols are identical to Fig. 1101 



are approximated as ^' = 127ri?^(_Bz), assuming a surface 
magnetic field Bs « 3(_Bz). The magnetic fluxes are plotted 
in Fig. [Inversus Lx together with the Pevtsov-relation and 
a sample of cTTS from Johns-Krull p007p . The tendency 
of the cTTS to show lower than expected X-ray luminosi- 
ties might be related to the effects of the disk, such as e.g. 
reduced coronal heating of mass-loaded magnetic field lines 
(Preibisch et al. I2005P or reduced height of coronal loops 
(Jardine et al. dOOSl)- Apart from HD 144668 (the lowest 
filled circle in Fig. [T^ , the Herbig Ae stars cluster around 
the line which follows the power law relation derived by 
Pevtsov. The lowest symbol corresponds to the debris disk 
star HD 172555. 

Clearly, the presence of the X-ray emission in Herbig 
Ae stars still remains to be explained and more spatially 
resolved observations are necessary, especially for stars with 
detected magnetic fields. 

5.3. Rotation 

Any kind of correlation between the magnetic field strength 
and the rotation period is of particular importance in view 
of the unknown origin of magnetic fields in Herbig Ae stars. 
For example, a trend of the magnetic field strength being 
lower in more slowly rotating stars would be consistent with 
the usual prediction for the dynamo theory for stars of a 
given mass (Mestel I1975|) . For no Herbig Ae star or debris 
disk star is the rotation period known. For stars with known 
disk inclinations and v sin i values we can estimate v^q val- 
ues, and from the knowledge of Veq and stellar radius we 
can estimate the rotation period. The v sin i values collected 
from the literature are listed in Col. 4 of Table [7] followed 
by values in Col. 5. Also most of the radii were col- 
lected from the literature. Only for three Herbig Ae stars, 
HD 101412, HD 135344B, and HD 179218, and two debris 
disk stars, HD 109573 and HD 181327, have the radii not 



Fig. 14. X-ray luminosity versus magnetic flux for our sam- 
ple of Herbig Ae stars, and the cTTS sample from Johns- 
Krull (|2007p . compared to the power-law relation derived 
for the Sun and active field stars by Pevtsov et al. ()2003p . 
cTTS arc represented by crosses. All other symbols are 
identical to Fig. [TOl 



been estimated in the past. For these stars, we estimated 
the radii using the Stefan-Boltzmann law. The informa- 
tion on bolometric luminosity and effective temperature of 
the three Herbig stars was found in van der Plas et al. 
(|2008p . while for HD 109573 we used the values presented 
by Debes et al. pTOg]) and for HD 18 1327 we used the val- 
ues presented by Chen et al. (|2006bp . Whenever the liter- 
ature sources permitted, the uncertainties associated with 
these parameters were included. Using Vcq and radii listed in 
Cols. 5 and 6 of Table[7l respectively, we estimated the rota- 
tion periods, which are fisted in Col. 7 of Table [T] Magnetic 
fields are detected in stars with a large range of rotation ve- 
locities, from 6 km s^^ up to 300 km s~^. No obvious trend 
of the strength of the longitudinal magnetic field with the 
rotation velocity or rotation period is identifiable in the sta- 
tistical tests. The star HD 101412 with the strongest lon- 
gitudinal magnetic field = 454±42 G) clearly stands 
out in both distributions as it is the most slowly rotating 
Herbig Ae star, with the longest rotation period of more 
than 17 days. It is possible that the slow rotation is caused 
by magnetic braking. 

5.4. Age 

There are clear indications for a trend towards stronger 
magnetic fields in younger Herbig Ae stars (Fig. [T5)) . con- 
firmed by statistical tests. The observations of magnetic 
fields of Herbig Ae stars, their strength and geometry are 
of a particular importance to understand the origin of mag- 
netic fields in Ap stars. It has been frequently mentioned 
in the literature that magnetic Herbig Ae stars are po- 
tential progenitors of the magnetic Ap stars (e.g., Stepien 
& Landstreet ?Mm Catala UDDl Wade et al. On 
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Fig. 15. The strength of the longitudinal magnetic field as 
a function of age. The symbols are identical to Fig. [TOl 

the other hand, from Fig. [T5] it is obvious that stronger 
magnetic fields appear in very young Herbig Ae stars, and 
the magnetic fields become very weak at the end of their 
PMS life. These results clearly confirm the conclusions of 
Hubrig ct al. (pOiM [MI51 I2007ap that magnetic fields in 
stars with masses less than 3 Mq are rarely found close to 
the ZAMS and that kG magnetic fields appear in A stars 
already evolved from the ZAMS. In contrast, magnetic Bp 
stars with masses M>3Mq seem to be concentrated closer 
to the ZAMS. 

5.5. Disk properties 

In the following, we compare the magnetic field and accre- 
tion rate to some relevant disk properties. 

5.5.1. Disk inclination 

Although an expanding sample of Herbig Ae stars have sys- 
tem inclination data from CO and coronographic imaging 
surveys, for a major part of the studied Herbig Ae stars 
the orientation of the disk still has to be constrained us- 
ing line profiles as a proxy for the inclination, as described 
in Sect. |4l The emission profile shapes of Mgii lines from 
UV studies or the Ha line profiles are frequently used. We 
present the inclination angles collected from different pa- 
pers in the second column of Table [T] Our own estimates 
of profile types from the study of Mgii spectral line pro- 
files in UV spectra are listed in the comments with running 
indices close to the HD numbers. In the third column we 
present the adopted inclinations. Using inclination angles 
in combination with usini data, wc derive v^q (Col. 4 of 
Tabic [21). 

However, no correlation between mass-accretion rate 
and disk inclination is detected. We also do not find any 
trend of the measured longitudinal magnetic field with the 
inclination angle fFig. [TO)) . 



Fig. 16. The strength of the longitudinal magnetic field 
versus disk inclination angles. Symbols are identical to 
Fig.Iini 



The disk inclinations for our magnetic Herbig Ac stars 
sample the whole possible range from close to face-on to 
close to edge-on. The strongest magnetic fields in our sam- 
ple are observed for a flared disk (HD 97048), for an edge-on 
disk (HD 101412) and for close to pole-on disks (HD 150193 
and HD 190073), and systems with an envelope rather than 
a disk (HD 85567 and HD 176386). These results show that 
there is no preferred disk orientation for the detection of a 
magnetic field. 

5.5.2. Disk geometry 

We find that magnetic fields are detected in systems both 
with and without strong PAHs, but appear to be more fre- 
quent in the first, flared case. This result could probably 
indicate an age dependence of the presence of magnetic 
fields. Acke et al. (|2004p studied Herbig Ae/Be stars with 
different IR SEDs with the conclusion that flaring disks 
probably evolve into self-shadowed disks. The existence of 
strong PAH emission may point to a flared disk, since the 
PAH grains reside further outside in the disk, whereas they 
would be destroyed near the hot inner disk rim. On the 
other hand, Keller et al. (|2008p conclude that most Herbig 
Ae/Be stars have PAH emission at some level, and there 
is no correlation of PAH emission with the disk geome- 
try. Since PAHs originate from the inner 200-300 AU of 
the disk, they are actually not a good measure of magnetic 
properties of the Herbig stars. Silicates usually emerge from 
~1AU and beyond. Contrary to the silicates, the PAHs 
need to be excited by direct stellar irradiation, which may 
not work well in a flat disk. In addition, the IR slope of the 
SED can further indicate whether the disk is flat or flared 
(Meeus et al. l20Qip . We note that magnetic fields arc also 
detected in HD 176386. 
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Table 7. Rotation periods of stars with known inclinations. 
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MecLSurements 
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L J 


[km s"^] 


[km s"^] 






XJT\ OCQQI 

tiu yooor 


/in 


/in 

4U— 00 


OU 


Dl — 10 


1 rfrf 
o.i 


on OR 


HD 97048 


30—40" 


30-40 


140'' 


218-280 


2.5'='= 


0.5—0.6 




42.8*" 












HD 100453 


>55^ 


>55 


39" 


<48 


1.7^^' 


>1.8 


HD 100546 


42'' 


47±5 


65* 


82-97 


1.7^^' 


0.9-1.0 




51"* 












HD 101412 


>55^ 


>55 


5" 


<6.1 


2 19S 


>17.4 


HD 135344B 


IV 


11 


69" 


362 


2.49s 


0.3 




14/ 












HD 139614 


<40'' 


<40 


15^^ 


>23 


1.6'*'* 


<3.5 


HD 144432 


48±10'= 


48±10 




83-114 


3.0-^ ^' 


1.3-1.8 


HD 144668 


58*^ 


58 


100™ 


118 


3.9'='= 


1.7 


HD 1 50193 


~38' 


30^ 


100' 


200 


2.1'='= 


0.5 






UO 1 u 


IS ^+ P 
lo.ozn _L 


±0.0 — z ± . 


9 fi'^'^ 


(S 1 —7 1 


HD 158643 


~90'"' 


90 


267±5" 


267±5 


5.3'"'' 


1.0 








228"" 


228 






HD 163296 


OJ-_9 


51 


130™ 


167 


2.8'='= 


0.8 


nu royr^z 


1 9 1 0^ 


1 1 


1 0^ 

DOzhz 


0/i'7 Qcr^i 
Z4/ — OOD 


l.D 


no HQ 


VVSer 


72±5''' 


72 


142^ 


149 


2.4'='= 


0.8 








229^= 


241 






HD 179218 


40±10'= 


40±10 


60'"" 


78-120 


4.8»» 


2.0-3.1 


HD 190073 


<40'* 


<40 




>18.7 


3.3^^ 


<8.9 


HD 9672 


90" 


90 


196^' 


196 


1.7" 


0.4 


HD 39060 


87" 


87 


ISO*" 


130 


1.8^'^' 


0.7 


HD 109573 


73° 


73 


152^" 


159 


1.8»s 


0.6 


HD 181327 


31.7^' 


31.7 


16'='= 


30 


1293 


2.0 



Notes: Inclination angles derived for resolved disk detections from the literature and from the study of Mgll spectral line profiles 
in UV spectra or the study of Ha line profiles in optical spectra and the inclination values we finally adopt for further use are 
listed in Cols. 2 and 3. v sin i values and radii presented in Cols. 4 and 6, respectively, are for the most part gathered from the 
literature. In Col. 7 we present the equatorial velocity values which are obtained from the inclination angles and the v sin i values. 
Vcq and radii are used to estimate the rotation periods listed in the last column. Radii marked as this study were calculated using 
Tcff and bolometric luminosity, both taken from the literature. 

References: "Doering et al. (pHoT)) . ''Doucet et al. p007| . '=Ardila et al. ([2007]! . ''Augereau et al. (pOOTj) . '=Dent et al. lf2005|l . 
^Pontoppidan et al. p008|) . »Meeus et al. lT998l ''Preibisch et al. p006|l . Tukagawa et al. p003)l . •'Alencar et al. (p003)) . 
^Slettebak (fT982|) . 'Wassell et al. (f2006ll . "Hug hes et al. (f2008|l . "Heap et al. (pOOOll . "Schneider et al . ([T999)l . ^'Schneider et al. 
((2006)) . ''Grady et al. (fT996ll . "Van den Ancker (fT998)), "Ack e & Waelkens ((2004)), *Hamidouche et al (pOOS)! . "this st udy (s ee 
Subsect.mTJ, "Dunkin et al. ffTWf} . "Hubrig et al. (|2007b|) . ^Royer et al. (gDDT}, ^Vieira et al. (^005)) , ^Mora et al. (pUDTj) . 
""Bernacca & Perinotto (fT970ll . ''''Pogodin et al. (|2005)l . '='=de la Reza & Pinzon (f2004)l . ''''Blondel & Tjin ([2006ll . '='=Hmenbrand et 
al. (fT992|) , ■''■''Wade et al. (pOOT)) . ^^this study (see table caption), ''''TatuUi et al. ((2008)) . "Hughes et al. ((2008)) . and ^Mi Folco et 
al. 1^^. 

Comments on the inclination angles deduced in this study, mainly from lUE (LWR and LWP cameras) through analysis of line 
profile types. Profile shapes are classified according to Beals (|195ip . See also the discussion in Subsect. 14.11 ^Mgll alternates 
between P Cygni type I (red emission, with blue-shifted absorption, lUE LWP 30772) and type HI (double emission, lUE LWP 
13082), ^the Ha profile shown by Meeus et al. ()2002p shows characteristic double-peaked Ho? emission, ^see the discussion in 
Subsect. SH *Mgn has a type I P Cygni profile (lUE LWP 25762, LWR 11977, LWR 08996), ^Mgii has a type I profile (lUE 
LWP 13083). 



5.5.3. Hot disk gas 

For many of our targets with a detected magnetic field, 
reports on the detection of excited, hot inner disk gas exist, 
which may be attributed to a gaseous accretion ring close 
to the star (see especially Tatulli et al. I2007( Isella et al. 

mm . 

Activity in the inner disk is often accompanied by an 
increased K-band excess. To prove whether there is a direct 
trend, we correlated the K-band excess, derived as the dif- 
ference of the measured and expected K-band flux of each 
star, with the strength of its magnetic field, but we found 



no correlation. We compared intrinsic V-K colours, as given 
e.g. by Koornneef ((1983p . to the difference V-K in the mea- 
sured photometry of both passbands, as given in SIMBAD. 
Due to a missing luminosity class wc assumed a class V 
for PDS 2, HD 97048, VV Ser and HD 179218. We obtained 
excesses between 0.0 and 3.9 mag, with an uncertainty of 
0.1 mag. The lack of correlation may be explained by the 
fact that without spectral resolution it is not possible to 
determine if the /•C-band excess is due to gas lines or to a 
hot inner dust continuum. A hot inner dust continuum will 
as well produce an increased broad-band K excess, possibly 
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without being related to the magnetic field. While there is a 
clear trend that most of the magnetic Herbig Ac stars have 
reports on hot, inner gas, a selection effect cannot be ruled 
out, as many of the non-magnetic Herbig Ae stars may just 
not have been measured in a similar way. 

5.6. Binarity 

A majority of the Herbig Ae/Be stars are binary or multi- 
ple systems. The existence of a magnetic field, however, 
does not seem to be related to binarity. We note how- 
ever, that the presence of a close companion contributing to 
the observed spectropolarimetric spectra can cause a non- 
detection of the magnetic field due to blending of spectral 
lines of the primary. The contribution of the secondary com- 
ponent can be disentangled only by means of high resolution 
spectropolarimetric observations, but not with low resolu- 
tion FORS 1 spectropolarimetry. In the present study, some 
targets, e.g. HD 97048 or HD 100546, show a magnetic field, 
but have no known companions, while other binary Herbig 
Ae stars possess no magnetic field. 

6. Summary 

In the course of our spectropolarimetric study of Herbig 
Ae/Be stars, we detected magnetic fields in six stars for 
the first time, PDS2, HD 97048, HD 100546, HD 135344B, 
HD 150193, and HD 176386. The presence of a magnetic 
field was confirmed in the stars HD 101412, HD 144668, and 
HD 190073. 

We have for the first time examined the relation of the 
measured field strengths to various parameters that char- 
acterize the star-disk system. 

Among the most important relations for the interpreta- 
tion of the fields is the one between magnetic field strength 
and accretion rate. We do not find a clear trend between 
these two parameters in our sample of Herbig Ae stars but 
the measured field strengths are compatible, in order of 
magnitude, with the values expected from magnetospheric 
accretion scenarios for a dipole field (tens to a few hundreds 
of Gauss). This contrasts with the situation for cTTS. For 
typical cTTS parameters, accretion models predict fields 
that range between ^^200-20000. However, the observed 
mean magnetic field strengths of cTTS are not correlated 
with the predictions (Johns-KruU [MI7|) . For most cTTS, 
the observed fields are larger than the expected values, pos- 
sibly indicating that magnetic field pressure dominates gas 
pressure in these systems. In addition, the dipole approxi- 
mation is known not to be valid for the case of cTTS that 
have complex field geometries (e.g. Gregory et al. 12008^ . 
while our results suggest that it may be a reasonable de- 
scription of Herbig Ac stars. 

We find that stronger magnetic fields tend to be found 
in younger Herbig stars. The magnetic fields become very 
weak or completely disappear in stars when they arrive on 
the ZAMS. Similarly, strong X-ray sources arc only found 
at the youngest ages, in qualitative agreement with the 
predictions of a shear dynamo that decays within a few 
Myrs as the rotational energy of the star decreases (Tout 
& Pringle fT995p . It is premature, however, to claim a direct 
connection between magnetic field and X-ray luminosity. 
The Herbig Ae stars seem to follow the power-law between 
magnetic flux and X-ray luminosity established for the Sun 
and main-sequence active stars. 



We do not find any trend between the presence of a 
magnetic field and disk inclination angles. The membership 
in binary or multiple systems does not seem to have any 
impact on the presence of a magnetic field, whereas there 
is a hint that the appearance of magnetic fields is more 
frequent in Herbig stars with flared disks and hot, inner 
gas. Since flared disks are the least evolved, this is possibly 
another indication for the decay of magnetic flelds with 
increasing age. However, no trend of the strength of the 
magnetic field with rotation velocity and rotation period 
was detected in our study. 

While considerable progress has been made with respect 
to the presence of magnetic fields in Herbig Ae stars, a num- 
ber of questions remain open. The most important question 
is related to the origin of the magnetic fields in these stars. 
Although our results provide new clues, the observational 
results presented in this work are still inconclusive as to 
the magnetic field origin. Tout & Pringle (|1995p proposed 
a non-solar dynamo that could operate in rapidly rotating 
A-type stars based on rotational shear energy. Their model 
predicts that the coronal activity at the observed rates of 
log Lx can be sustained for a period of the order of 10^ yr. 
Other possible mechanisms causing magnetic activity in- 
volve fossil magnetic fields or magnetically confined wind 
shocks (e.g. Babel & Montmerle ll997p . A more comprehen- 
sive survey of the presence of magnetic fields and a detailed 
study of the magnetic field topology in a Herbig star sample 
of increased size will provide important additional informa- 
tion to test the predictions of different theories. About half 
of the stars with magnetic field detections possess longitu- 
dinal magnetic fields larger than 100 G. These stars are the 
best candidates for future spectropolarimetric studies to 
analyze the behaviour of their magnetic fields over the ro- 
tational periods to disclose the magnetic topology on their 
surfaces and to study the complex interaction between the 
stellar magnetic field, the disk and the stellar wind. 
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